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Abstract. We propose a scheme in which the quantum coherence of a 
nanomechanical resonator can be probed using a superconducting qubit. We consider 
a mechanical resonator coupled capacitively to a Cooper-pair box and assume that 
the superconducting qubit is tuned to the degeneracy point so that its coherence 
time is maximised and the electro-mechanical coupling can be approximated by a 
dispersive Hamiltonian. When the qubit is prepared in a superposition of states 
this drives the mechanical resonator progressively into a superposition which in turn 
leads to apparent decoherence of the qubit. Applying a suitable control pulse to 
the qubit allows its population to be inverted resulting in a reversal of the resonator 
dynamics. However, the resonator's interactions with its environment mean that the 
dynamics is not completely reversible. We show that this irreversibility is largely due 
to the decoherence of the mechanical resonator and can be inferred from appropriate 
measurements on the qubit alone. Using estimates for the parameters involved based 
on a specific realization of the system we show that it should be possible to carry out 
this scheme with existing device technology. 
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1. Introduction 

One way of exploring tlie quantum colierence properties of a nanomechanical resonator 
is to couple it to a qubit formed by a solid-state two- level system (TLS). Coupling to 
an isolated harmonic oscillator can initially cause an apparent loss of phase coherence 
in the qubit if the oscillator is driven into a superposition of orthogonal states, but 
signatures of the overall coherence of the full system (i.e. oscillator and TLS together) 
can be found in the subsequent dynamics of the TLS. However, if instead the qubit is 
coupled to a harmonic oscillator which is in turn coupled to a bath, then the effective 
dynamics of the TLS and oscillator will now be different and the loss of the oscillator's 
coherence due to the bath will be manifest in the dynamics of the TLS [H [21 [3]. 

From a theoretical point of view it is relatively straightforward to devise simple 
schemes based on these principles to probe the rate at which the environment causes 
decoherence of an oscillator [3] . Indeed, exactly this kind of approach has been used 
very successfully in the field of cavity quantum electrodynamics (cQED) to probe the 
quantum coherence of a mode of the electromagnetic field by examining its influence 
on effectively two- level atoms (5} El [7] . Similar experiments have also been carried out 
successfully on trapped ions with the internal electronic state of the ion playing the role 
of the two-level system and the ion's motional state the oscillator [8], [9], ITO] . 

The development of relatively large and well controlled quantum coherent TLSs 
in the solid-state, such as superconducting circuits designed to act as qubits, seems to 
offer a way to perform analogous experiments with nanomechanical resonators [H [2]. 
Furthermore, recent experiments have demonstrated that it is possible to recreate many 
of the features of traditional optical cQED in the solid state using a superconducting 
qubit coupled to a superconducting resonator pTl [12] . Since nanomechanical resonators 
are typically a few microns in length and contain macroscopic numbers of atoms, 
producing a quantum superposition of spatially separated states in such systems and 
monitoring its progressive loss of coherence (due to interactions with its environment) 
would represent an important increase in size of the system involved compared to 
superpositions of ions and light [TSj [T3]. However, performing quantum coherent 
experiments using a nanomechanical resonator is likely to be more difficult than with 
a superconducting one as nanomechanical resonators are generally much lower in 
frequency. 

In order to understand the practical difficulties entailed in using a superconducting 
qubit to probe the decoherence of a nanomechanical resonator we briefly review 
the apparent constraints which any such scheme must satisfy. First of all, the 
superconducting qubit must remain sufficiently coherent that the influence of the 
mechanical resonator's environment can be clearly discerned in its dynamics. Secondly, 
it will be desirable to couple the TLS and resonator as strongly as possible since the 
signal (s) of coherence and/or decoherence in the mechanical resonator measurable in 
the TLS will become clearer the larger the coupling between the TLS and resonator. 
Finally, unless impressive cooling of the resonator can be achieved, the experiments will 
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always have to contend with the competing effect of phase smearing arising from the 
range of oscillator states (and their associated phases) in the thermal ensemble of the 
oscillator. Again the relatively low frequencies of mechanical resonators make this more 
of a problem than it would be in the superconducting case. Note that in practice there 
is no simple way of designing a double clamped beam resonator to optimize all of these 
constraints at onc^. For example, in most realizations the TLS-resonator coupling will 
increase as the resonator is made larger, but enlargement of the resonator will inevitably 
reduce its fundamental frequency. 

In this paper we describe how a dispersive interaction between a superconducting 
qubit and a nanomechanical resonator can be used to produce superpositions of the 
resonator state and how the coherence of this superposition can then be probed by 
measuring the state of the TLS. We identify the relatively short coherence times of the 
superconducting qubit as the most serious constraint on these types of scheme and hence 
assume that the TLS is tuned to operate at a point where its coherence is maximized; 
it is this choice of operating point which leads to a dispersive coupling between the 
TLS and the resonator. Although the dispersive interaction is relatively weak, we find 
that the effect of the TLS on the resonator can be amplified by preparing the latter 
in a state with large amplitude. We explore in detail the quantum dynamics of the 
resonator and TLS including the effects of the inevitably mixed initial state of the 
resonator and the interaction with the environment. In assessing the extent to which 
the schemes we propose can be carried out in practice we make use of the analysis in 
the companion reference [16], which considers how sufficiently strong coupling between 
a nanomechanical resonator and a superconducting qubit can be best achieved without 
degrading the coherence of the qubit. 

Our work builds on and extends previous studies of similar systems [U |2[ [T71 [3l [18]. 
In particular, we believe that the scheme outlined here represents an important 
improvement on that proposed by us in Ref. [1] in a number of respects. Most 
importantly, the scheme we propose here is more likely to be practicable as it is designed 
to be performed at the degeneracy point of the qubit where it remains coherent for 
at least an order of magnitude longer [T9] than the operating point considered in [1]. 
Furthermore, the scheme is much more flexible in the sense that it would be possible 
to vary several of the important parameters systematically (such as the phase space 
separation of the resonator states involved and the duration of the superposition). This 
would be an important advantage in interpreting the results of this type of experiment, 
since the nature of the mechanical resonator's environment is not well understood [20] 
and in a sense the purpose of the experiment we propose would be to provide empirical 
information about it. 

This paper is organised as follows. In section[2]we introduce the generic Hamiltonian 

I One interesting way to avoid the problems posed by the relatively low frequency of flexural-mode 
mechanical resonators is to use a different type of mechanical system, such as dilational disk resonators 
which can have frequencies well beyond 1 GHz [15 . Here, however, we will confine our attention to the 
conventional mechanical resonators formed by doubly clamped beams. 
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for the superconducting qubit-resonator systems which we will work with here. We 
discuss the practical constraints which dictate our choice of operating regime and 
introduce the effective (dispersive) Hamiltonian which is valid when the mechanical 
resonator is much slower than the superconducting circuit. Next in section [3] we describe 
how the dispersive interaction can create states which involve superpositions of spatially 
separated states of the mechanical resonator so that measurements on the TLS alone 
show an apparent initial loss of coherence. We show that the TLS coherence can be 
recovered (recoherence) in a controlled way using a particular choice of control pulses. 
In section H] we calculate how the presence of the environment of both the mechanical 
resonator and the qubit itself affects recoherence. We then consider the values of 
the various parameters which are likely to be practicable in present or near future 
experiments in section [H Then in section [6l we present calculations of the behaviour of 
the recoherences for a range of practicable parameter values. Finally, section [7] contains 
a discussion of our results and our conclusions. The Appendix contains further details 
on some of our calculations. 

2. Resonator-TLS Effective Hamiltonian 

2.1. Operating Regime 

The generic Hamiltonian for the superconducting TLS and mechanical resonator which 
we consider here is, 

H = ^cr^ + Acr^ + \{a^ + a)a^ + hw ^a^a + (1) 

where the qubit energy scales eo and A depend on the details of the specific 
superconducting system considered, u is the resonator frequency, A is the strength of 
the resonator-TLS coupling, and the operators az(x) act on the TLS while a'-^-' act on the 
resonator. The mechanical resonator is assumed to be the fundamental flexural mode of 
a suspended doubly-clamped beam. The coupling between the TLS and the mechanical 
resonator is implicitly assumed to be weak in the sense that only linear (in the resonator 
position coordinate) coupling needs to be considered. The TLS states are defined as |1) 
and |0) so that e.g. = |1)(0| + |0)(1|. 

This Hamiltonian is derived in Ref. [T6] for the specific cases of either a Cooper 
pair box (CPB) with an island which is suspended to form the mechanical resonator, or 
a flux qubit (again with suspended segment that forms the mechanical resonator). In 
each case the qubit is also assumed to be fabricated close to the centre electrode of a 
superconducting microwave coplanar waveguide (CPW) resonator. The CPW resonator 
provides a way to both measure and manipulate the qubit state, as well as a means 
to drive the mechanical resonator into an initial state which has a large amplitude. 
Although both the mechanical resonator and the qubit are also coupled to the CPW 
resonator in this system as just stated, we will not include the latter explicitly in this 
article, assuming that it is unpopulated (i.e. it is at or close to the vacuum state), except 
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initially when used to drive the mechanical resonator [TH] and for the short periods when 
it is used to manipulate the state of the qubit. Population of the CPW resonator is 
required when measuring the qubit state, but at this stage disruption to the mechanical 
system and dephasing of the qubit are unimportant so long as the measurement can still 
be performed with high fidelity, which we will assume is the case. For a CPB the states 
|0), |1) correspond to different charge states and the coupling between the TLS and the 
resonator is capacitive. In contrast, when the Hamiltonian [equation (Q] is realized with 
a flux qubit, the coupling between the TLS and the mechanical resonator is inductive 
and the relevant qubit states are of current circulating in opposite directions [2| [T6l [2T] . 

The best coherence times for both superconducting charge and flux based qubits 
are achieved at the degeneracy point where eo = 0. Away from the degeneracy point, 
experiments [19] have demonstrated that the coherence times of superconducting qubits 
decrease by orders of magnitude. We regard the coherence of the superconducting TLS 
as the primary constraint and hence we choose to operate the TLS at its degeneracy 
point when probing the resonator's coherence. Another important constraint arising 
from the use of the superconducting TLS is the need to avoid thermal mixing of the two 
states involved. In practice, for experiments performed at temperatures of order 20 mK 
this means that we will require energy separations between the two states (i.e. 2A when 
working at the degeneracy point) in the superconducting TLS that are much larger than 
the thermal energy scale. Experiments [22] using a CPB (coupled to a superconducting 
CPW resonator for state control and read-out) achieved coherence (T2) times of up to 
0.5 /is (operating at the degeneracy point) and relaxation times (Ti) of about 7 fis with 
z/q = 2A/h ~ 5 GHz and we take these as indicative of the current practical limitations. 
Note that even longer coherence times of up to 2 fis have been reported for experiments 
on CPBs using an echo technique in which the TLS state is inverted mid-way through 
the experiment [23] 

In terms of the mechanical resonator, we will consider beam structures which are 
fabricated by under-etching a bulk substrate or metallic film. The fundamental (fiexural) 
mode frequencies of such devices can in practice be as high as 1 GHz [21], but because the 
electro-mechanical coupling. A, (in equation [1]) for such modes increases with the length 
of the beam, /, (for both capacitive and inductive couplings [16]) whereas the frequency 
clearly decreases with increasing /, it is clear that high frequencies can only be achieved 
at the expense of very weak couplings [T6l [H [2] . Nanomechanical resonators have 
already been fabricated in close proximity to superconducting structures [25], but with 
mechanical frequencies ~ 20 MHz. We therefore assume that the mechanical frequency 
will be much lower than the energy scale of the qubit, i.e. A ^ huj. Having made this 
assumption of a wide separation of time-scales for the mechanical and superconducting 
elements we can proceed to derive a simpler effective Hamiltonian which is valid in this 
regime |§| 

§ Note that the CPW resonator used to manipulate the CPB state will have a frequency which is close 
to Va and hence is also very far from the mechanical frequency. 
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2.2. Adiahatic limit 

When eo, is tuned to zero it is convenient to rewrite the Hamiltonian in terms of a new 
basis for the qubit. Defining new basis states, 

l±) = ^(±|0) + |l)), (2) 

we can write the Hamiltonian ([1]) for e = as 

ifdcg = + A(a^ + ayUx + hw (^^a + (3) 

where the new spin operators are defined in terms of the new basis states. 

We proceed by exploiting the separation in time-scales to make an adiabatic 
approximation [21 EZl EB|- Since the mechanical resonator will generally be in a 
Gaussian state of large amplitude, in what follows it is reasonable to take a semi- 
classical approach [26] . We initially assume that the mechanical resonator is at a fixed 
position X and use this to calculate the eigenvalues of the TLS, these are then used to 
calculate an effective Hamiltonian for the oscillator. With the resonator at position x 
the Hamiltonian of the TLS is (using equation [3]): 

Htls = ^'^z + K^/XzpP^,, (4) 

where x^p = {h/2mujY^'^. The eigenvalues of the TLS are now e± = i-y/A^ -|- (Xx/x^p)^- 
For sufficiently weak coupling i.e. [Xx / {x zpA)]"^ <^ 1, we can expand the eigenvalues to 
lowest order 

e^ = ±A{l + ^{-^{ ]. (5) 




The evolution of the mechanical system over time then causes a weak position dependent 
(and hence ultimately time dependent) perturbation to the eigenstates of the TLS. 
In the adiabatic approximation the wide separation of time-scales and weak coupling 
mean that the TLS evolves smoothly within each eigenstate with its dynamics arising 
from changes in time of the eigenstates themselves, rather than any transitions between 
different eigenstates. 

Labeling the instantaneous eigenstates of equation (jl]) as |+) and |— ), we can write 
down the effective Hamiltonian felt by the oscillator for the TLS confined to one of its 
eigenstates as, 

Hi = hw (^ata + ± (^A + ^(at + af^ (6) 

Therefore within the framework of the adiabatic approximation we can write down the 
following model Hamiltonian for the system, 

H = A(l + ^{a^ + a)Aaz + nw(a^a + l]. (7) 



We have dropped the distinction between perturbed and unperturbed eigenstates as it 
does not play a role in what follows. 
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Assuming that the couphng term is a weak perturbation (i.e. assuming A^/2A -C 
hw), we can also make the rotating wave approximation in which the terms and (a^)^ 
are dropped [2]. The final result of these approximations is the following dispersive 
Hamiltonian for the TLS-resonator system, 



where ui = A^/ (hA). A key feature of this Hamiltonian is that the perturbation of the 
oscillator commutes with the unperturbed Hamiltonian (i.e. it is a QND Hamiltonian). 
This feature is exploited in schemes to measure the number state of a resonator using 
a superconducting qubit [SOj [30l [31] and it also plays an important role in what follows 
here. Note that this Hamiltonian can also be obtained from equation ([3]) via a range of 
other approaches [H [HI [32l |30l [33] . 

3. Coherent oscillations and recoherences: simple description 

The dispersive Hamiltonian shifts the mechanical frequency in a way that depends on 
the state of the TLS. This interaction can be used to probe the quantum coherence of 
the mechanical resonator. The idea is to perform a Ramsey interference [7] experiment 
in which the qubit is prepared in a superposition of its eigenstates states using a control 
pulse, this superposition is then allowed to interact with the resonator for a time t 
before a second pulse is applied to the qubit and then a measurement of its state is 
performed. For an isolated TLS the probability of finding the system in one or other 
of its eigenstates at the end of the experiment will oscillate between zero and unity as 
a function of the time between the two control pulses. When the mechanical resonator 
is present the interaction with the superposition of TLS states leads to an overall 
superposition of states involving spatially separated mechanical states. For a sufficiently 
strong interaction, the separation of the resonator states coupled to the qubit states leads 
to a strong suppression of the oscillations in the final qubit state measurements. The 
coherence of the resonator can be inferred by inverting the state of the TLS midway 
between the two original control pulses. The scheme is illustrated schematically in figure 
[1] In the absence of the resonator's environment, such an inversion should lead to a 
reversal of its dynamics and hence the recovery of the oscillations in the final TLS state 
measurement J6\. Very similar schemes have been demonstrated in optical systems [7]. 

In what follows, we will assume that it is possible to measure the state of the TLS 
within the |±) basis and to rotate its state by applying transformations of the form 
exp{—i6ax/2) with a parameter 6 that can be controlled to a high degree of precision. 
These requirements are readily met in the system of a charge or fiux qubit (with 
suspended segment forming a mechanical resonator) coupled to a superconducting CPW 
resonator described in Ref. |T6], which we consider here. The TLS state is determined 
by measuring the transmission of an off-resonant pulse applied to the CPW resonator, 
while rotation of the TLS is performed by applying almost-resonant pulses to the CPW 
resonator and making use of the resulting Rabi oscillations [32l 122] . 
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We will begin by considering the simple though unrealistic case of an isolated 
resonator which is initially prepared in a coherent state, |ao). The effects of the 
environment on the evolution of the coupled TLS-resonator device and the types of 
initial (mixtures of) states of the oscillator that can be prepared in practice are addressed 
in later sections. We assume that the TLS system is in its ground state |— ), hence the 
total initial state is |— ) (S) \ao)- Application of an appropriate control pulse to rotate 
the state of the TLS by = 7r/2 produces a superposition of TLS states. Since the 
rotation of the TLS will in practice be very fast compared to the mechanical period we 
can neglect any evolution of the mechanical resonator during the pulse and hence write 
the total state of the system after the pulse as p(0) = \ip{0)){ilj{0)\ with 

IV'(0)) = ^(|-)-.|+))®|ao), (9) 

Starting with this initial state at t = 0, the dispersive interaction [equation (IHl)] leads 
to the following joint state after time t, 

\m) = 7f ® - ^1+) ® (10) 

where 

|a±(t)) = e^*'^*/'^e-'("±"^)("^"+^/2)t|^^^ (^^^ 

The resonator evolution in phase space during this period is illustrated in figure [T]d. The 
next step is to perform a second 7t/2 rotation on the TLS, leading to the state 

iV^W(t)) = l[\-)® W) - Mm -^\+)® W) + Mt)))] .(13) 

Finally, the state of the TLS is measured in the |±) basis. The probability of finding 
the TLS in state |+) for a period of evolution t between the two control pulses is 

ii+)(t) = Tr[|+)(+|p(t)] = i (1 + Re[{a4t)\a4m) . (14) 
The overlap is readily evaluated, 

(a_(t)|«+(t)) = e-\-0\'{l-e-'^-^^)^-2^At/H-^.,t_ (^5^ 

The final result for P\+){t) is thus, 

i^^^(t) = ^ + ^Re [e-l"ol'(i-^""^^')e-2*^*/^-*-^*] . (16) 

Note that this function depends only on the amplitude of the initial mechanical state, 
not its phase. 

The behaviour of P\+){t) is easy to understand. Without any coupling to the 
resonator the coherent oscillations in the TLS state mean that the probability oscillates 
over time between zero and unity with a period tr = h/{2Aji^ a key indicator of the 

II Note that in practice the pulses used to rotate the state of the TLS are chosen to be shghtly ofF- 
resonant. As a resuh a stroboscopic observation of the osciUations in P|_|_^ can be made which replaces 
the very fast oscillations at frequency 2A/h with much slower (and hence easier to observe) ones at the 
chosen de-tuning frequency [7l 1341 122j , we neglect this detail here as our primary interest is not in the 
frequency of the oscillations but in their amplitude. 
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Figure 1. Schematic illustration of the evolution of the mechanical resonator in phase 
space during the echo sequence. Initially (a) the resonator is prepared in a coherent 
state and the qubit is prepared in a superposition of states. The two qubit states 
couple to the resonator leading to different effective frequencies w ± wi so that in the 
frame rotating at the resonator frequency the two mechanical states start to pull apart 
(b). A TT pulse inverts the qubit state and hence interchanges the relative frequencies 
of the two resonator states (c). When the periods of evolution before and after the 
inversion of the qubit are the same the resonator will return to its initial state (d) in 
the absence of dissipation. 



quantum coherence of the TLS [22] . For sufficiently strong couphng, the resonator causes 
a relatively rapid reduction in the amplitude of the oscillations as a function of time 
leading to a period where = 0.5, implying that the resonator decoheres the 

TLS. However, because the resonator is a periodic system and is itself coherent, the 
oscillations in P\+)(t) reappear, giving rise to so-called recoherence, for t ~ Tr/cui [2]. 

Although recoherence does occur naturally after a time t ~ 7r/u;i it is preferable to 
use an approach where the time between coherences can be varied systematically. This 
is readily achieved using a spin echo technique to induce recoherence at a chosen time. 
This type of approach was used with great success in optical cQED experiments [6j as 
well as experiments on superconducting circuits [T2| [23]. 



3.1. Echo technique 

For the spin echo sequence we again start with the system in the state \ip{0)) (equation 
[9]) and allow it to evolve as before for a time ti so that. 



mti)) = ^{\-) 



(17) 



Next we apply a control pulse to the TLS which effectively inverts the populations of the 
two eigenstates [this corresponds to the unitary operation exp{—i9ax/2) with 9 = n]. 
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Thus, just after the pulse we have 

mt)) = ^ (^1+) + I-) ® |a+(ti))) . (18) 

We now allow the system to evolve for a further time t2, after which the resulting state 
of the system will be 

+ t^)) = ^ ® |a_+(t2,ti)) + «|+) ® |a+-(t2,ti))) . (19) 



V2 



where now, 



|a_+(t2, ti)) = e-*^(*i-*2)/^e-*(^-"^)("^"+^/2)t2g-^(..+-i)(ata+i/2)ti 1^^^ (20) 
|a+_(t2, ti)) = e^^(*i-*2)/ V*("+^i)('^^'^+^/2)*^e~*("~^i)(''^'*+^/2)*^ |ao)- (21) 

Note that the simplicity of this expression relies on the fact that the perturbed resonator 
Hamiltonian commutes with the unperturbed on^, thus we find 

(a_+(t2,ti)|a+_(t2,ti)) = e2*^(*i-*^)/^(ao|e*'"'^"'"'''/'^^*'"*'Vo)- (22) 

Carrying out a final rotation of the TLS state (with 6 = 7r/2) the final overall probability 
of finding it in state |+) is given by 

P|^^(t, + t2) = ^- ^Re [e-l"ol'(i-^"^^^*^-'^')e2*^(*^-*^)/^+^-^(*^^*^)] . (23) 

The probability P\+) is zero at t = ti + 12 when ti = t2, this is because at this instant 
the oscillator states associated with each of the qubit states are the same so that the 
effect of the pulses is simply to rotate the qubit through a total of 27r. To examine the 
apparent coherence of the qubit we can define the envelope of the oscillations in P\+), 

E[Pi^){h + h)] = I + ie-^^[l-P(i-e^-<*-'^')l (24) 

= 1 _|_ g|aoP{l-cos(wi[ti-t2])}j _ (25) 

The envelope of the oscillations is unity when t2 = ti signifying the recoherence of the 
qubit. Thus we can use the echo approach to induce recoherences in the qubit dynamics 
whenever we choose by tuning ti(= ^2)- We note that this particular approach also has 
the advantage that inverting the state of the TLS at t = ti can lead to an increase in 
the effective coherence of the TLS as measured at ^2 — ti as it eliminates dephasing 
effects arising from fluctuations in the TLS energy level spacings which occur between 
different experimental runs [23] . 



3.2. State separation and entanglement 

After evolving for a time t (and without any inversion of the TLS states), the two 
coherent states of the resonator that (together with the TLS states) form a superposition 

% Note, however, that with an appropriate choice of control pulse at i = ti the idea of an echo 
experiment is not limited to systems with dispersive Hamiltonians. This type of experiment has been 
performed for systems with a Jaynes-Cummings type interaction 6J. 
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are \a±(t)) [equation (IT^ ] and they have a separation in phase space given by 

S{t) = [((Xi)., - {X,)^S + {{X,}^^ - {X,)^_rY^' (26) 

= 2|ao| sincjit (27) 

where {Xi)a^ = {a+\Xi\a+) etc, and the phase space operators are defined by 
Xi = (a + at)/2 and X2 = (a — a^)/(2i). Because we are deahng with a pure state 
of the TLS and resonator we can also obtain the entanglement of the system, E{t), 
by calculating the von Neumann entropy of the reduced density matrix (of either the 
resonator or the TLS)EI. Evaluating this we find that it is entirely determined by the 
phase space separation of the resonator states [13j . 

E{t) = 1 - log, [(1 + xf'-^'>/\l - x)^'-^^/'] (28) 
where x(t) = exp(— S'(t)^/2). The entanglement E(t) rapidly saturates at its maximum 
value of unity as the separation S is increased: it reaches about 0.75 for S* = 1 and is 
already very close to unity for = 2. Note that the decay of the qubit oscillations also 
depends on S alone: we can rewrite equation (fT6ll as, 

P\+){t) = ^{1 + X{t) cos ^{t)) (29) 

where the (real) phase is defined as (f){t) = {2A/h + uji)t + |a;oP sin(2a;it). 

The aim of these simple calculations is just to show that the separation of the 
resonator states in phase space provides an important figure of merit for the kind of 
experiment we have in mind. The diameter of the 'uncertainty circle' [33] for a coherent 
state is 1/2 and so one basic (though somewhat arbitrary) criterion for producing a 
distinguishable superposition of resonator states is to require S{t) > 1. Although 
according to equation (127|) the largest separation is achieved when uit = 7r/2, in practice 
the limited coherence times available for the TLS mean that the evolution times will 
be such that Uit <^ 1 and hence we can approximate S{t) ~ 2|ao|^i^- If we use the 
spin echo approach then the maximum separation will be achieved at t = ti just before 
the TLS state is inverted and hence to achieve a meaningful superposition we would 
need to have 2|q;o|i^i^i ^ 1- The details of how a driven resonator state could be 
prepared in practice for the qubit-mechanical resonator system in which the mechanical 
component is formed by suspending part of the qubit circuit is considered in reference 
[16j ; we will make use of the results obtained there when considering what kind of initial 
mechanical state could be prepared in practice, but for now we point out the crucial 
role played by the magnitude of the initial coherent state, |ao|. The size of resonator 
state superposition produced depends through ao on the initial state of the oscillator. 
This provides us with a way of overcoming the weak TLS-resonator coupling and the 
wide separation of their time scales: by preparing the resonator in a state with large 
enough |q;o| we can overcome the very weak interaction with the qubit to nevertheless 
produce relatively large superpositions over the relatively short times during which the 
TLS remains coherent. On the other hand, unless we start with a state with non-zero 
ao then S will be zero throughout (and no entanglement will be produced). 

Note that the entanglement dynamics of this system was studied very recently for mixed states [35j 
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4. Role of the resonator's environment 

We now consider tlie effect wliich the resonator's environment has on the recoherences 
in the qubit. The interaction between the resonator and its surroundings is typicaUy 
modelled by including a bath of oscillators that are weakly coupled to the resonator. 
This approach is the one followed in quantum optics and although it is not clear to 
what extent it represents an actual nanomechanical resonator's environment [20], it can 
at least be justified in the idealised case where dissipation in the collective mechanical 
mode which forms the resonator is due only to coupling to the bulk phonon modes in 
the supports of the resonator [36] . In this simplified description the effects of the bath 
on the 'system' resonator can parameterized by a damping rate 7 for the resonator and 
a temperature T^, which can be expressed in terms of the average number of the quanta 
the resonator would have if it were in equilibrium with the bath, 

^= eW-l - 

For sufficiently high temperatures (fce^r ^ ^7) the master equation for the mechanical 
resonator and qubit including the dissipative effects of the resonator's environment can 
be modelled using the quantum optical damping kernel [37] . 

P = Cp = -'-[H,,p] + C,p (31) 

where 

£.P--.^|:^.{p.ri]- "'""'7^^" K[^.p]]. (32) 

Assuming that the oscillator damping is very weak (7 <^ 00) we can further simply the 
dissipative part of the master equation by using the rotating wave approximation, 

CdP = —-^ [a^ap + pa'^a — 2apa^) — 7n [a^ap + paa^ — apa^ — a'^ pa) . (33) 

We stress again that we use this damping kernel here to provide a simple illustrative 
estimate of the dissipative dynamics of the mechanical resonator. The true form of the 
mechanical damping kernel remains somewhat uncertain and one of the aims of the 
experiments we propose would be to obtain empirical information about it. 

The superconducting qubit is also subject to decoherence due to interactions with 
other degrees of freedom in the system apart from the mechanical resonator . The 
dissipative dynamics of such systems can be characterised by the relaxation times in the 
equations of motion for the diagonal and off-diagonal components of the TLS density 
operator. The decay of the excited state population of the TLS is described by Ti, while 
the decay of the TLS coherence is described by T2. In practice, Ti times have been 
typically an order of magnitude larger than T2 times [22j . Since we will only consider 
total evolution times t (before measurement) of the system that are shorter than T2, we 
therefore will always have t <^Ti and hence can neglect relaxation of the TLS in what 
follows. The master equation for the system [equation (l3T|l ] can be written in terms of 
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P++ 


= A 


+P 






-P 


p+- 


= A 


-P 



the components p^^ = ) etc, incorporating a finite T2 time as follows, 

= -iuj+[a^a, P++] + Cdp++ (34) 

= — iu;_[a^a,p ] + CdP (35) 

= - {i2A/h + iuJi+T2^) p+_ (36) 
- iuj[a^a, p+_] - zco'i{a'''a, p+_} + CdP+-, 

where uj± = uj 

We assume that immediately after the first control pulse is applied to the TLS the 
state of the system is given by 

p{o) = \m){m\®pi;:°\ m 

where \ip{0)) = (|— ) — i\+))/V2 and pj^°'' is a displaced thermal density operator [38] 
for the resonator defined by 

p£°) = D{ao)pthD\ao) (38) 

d^u = (39) 

Tim 

where D{a) = exp(aa^ — a* a) is the displacement operator and we have defined 
z/ = a + ao in the last line. The undisplaced thermal density operator is 





e 



-\a\''/m 



Pth= / da ^—\a){a\, (40) 



nm 



where m = (e^^^'^^ — l) ^ . We have chosen to specify a temperature Tj for the initial 
state of the mechanics resonator which can be different from that of the environment 
Tr- Simply driving the resonator (assuming a noiseless drive) would ideally lead to a 
displaced thermal state with Tj = T,.. However, it is interesting conceptually to consider 
the case where the mechanical resonator is somehow pre-cooled to a lower temperature 
than its surroundings Tj < T^. Alternatively a choice of Tj > provides a simple model 
for the case where there is no cooling and instead the drive adds noise to the resonator 
state. Although the initial resonator state will be prepared by driving, we assume that 
the drive is switched off before the first pulse is applied to the TLS. 

The evolution of the component equations (1M|1 - (136|1 can be calculated very 
conveniently using a phase space approach [391 HHl EHl [301 SI]- The method involves 
working with the Wigner transform of the components defined as 

W+.{x,p;t) = —J dy{x + y\p+4t)\x-y)e-'^'''y/'' (41) 

etc, which evolve according to the set of (uncoupled) partial differential equations 
obtained by transforming equations ( 1MI) -( 1361) . For our choice of an initial displaced 
thermal state, each of the initial Wigner functions is Gaussian and remains so during 
the evolution. This means that the relevant partial differential equations for the Wigner 
function components can be solved via a Gaussian ansatz. Details of the calculation 
(which follows the approach used in reference ^0]) are given in the Appendix. 
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Using the phase space approach, we readily obtain the following expression for 



where 



with 



e{t) = - {i^/2 + uJip)t -i\n 



1-M 



a, 



— t- 



0\_ (^^~i2ujif3t _ 



I _ ]\/[Q-2iuJi/3t 

1-M 
I _ ][fQ-2iuJi/3t 



(3 = ([I - i-f/2ujif - 2i-fn/uj-^Y^^ 

M= (2m+l)-/?-^7/2u;i 
(2m + l)+/3-i7/2cJi' 

Note that in the limit 7 — we recover the much simpler expression [2] 



i^^^)(t) = i + le-*/^^Re 



-r){t)\ao\^/{l+m'n{t))^-i{2l\/h+w^)t 

1 + mrjit) 



with ri{t) = 1 — e^ 



-2iLu-i t 



(42) 
(43) 



(44) 
(45) 

(46) 



4.1. Echo sequence 

We now consider the case where an additional vr pulse is applied to the system at 
time t = ti after the first 7r/2 pulse, and then the final 7r/2 pulse is applied at time 
t = tf = ti + 12- The evolution of the density matrix between the two tt/2 pulses can 
be written as [121 E] 

p(tj) = e^(*^-*^)7^e^*V(0), (47) 

where 

Up = e-*^^-/2pe*^^^/2. (48) 

In order to calculate P|+)(t/) we need the off-diagonal component of the density matrix 
given by, 

p+^tf) = e^+-*V-+(ii) = e^+-*Vi-(ii)- (49) 

This evolution can again be calculated using a phase space approach (see the Appendix 
for details). The resulting final probability for finding the TLS in state |+) takes a very 
similar form to before, 

although the expression for 6{tf) is rather complicated [it is given in full in equation 
flA.36p ]. It is important to note that even though the system is damped, the phase 
space separation between the components of the mechanical resonator's density matrix 
corresponding to the diagonal elements of the TLS still vanishes at t = ti +^2 for ^2 = ti- 
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The use of an echo technique allows us to filter out many of the effects that arise 
just because we start with a mixed state such as a decay in the oscillation amplitude 
due to averaging over the different phases of oscillation associated with the different 
resonator states in the initial mixture. The recoherence 'signal' measured at the echo 
time t = 2ti is the irreversibility of the system's dynamics [12] . What we are in effect 
measuring is the dynamics due to the resonator's damping kernel. There is no simple 
way of partitioning the dissipation into a contribution from the decoherence of spatially 
separated states and simple fluctuations in the resonator's energy during the experiment: 
both contribute to what is measured. An important consequence of this is that a perfect 
recoherence is not achieved for 7 7^ even if ao = 0. However, when relatively large 
phase space separations of the resonator state are achieved (S* > 1) and the experiment 
is performed on a time-scale which is very short compared to the energy relaxation time 
1/7, we can expect the decoherence of the superposition of mechanical states to be the 
dominant contribution to the irreversibility of the dynamics. 

5. Practical considerations 

We now turn to the question of what kinds of parameters might be achievable in practice 
and hence the prospects for using the approach we have been discussing to probe the 
quantum coherence of a nanomechanical resonator in the near future. A key quantity 
which we need to examine is the maximum phase space separation, S{ti), between 
resonator states that can be achieved at the mid-point of an echo experiment. As we 
have seen, a large initial amplitude for the resonator |ao| will enhance the phase space 
separation. However, for our theoretical approach to be valid we need to ensure not 
just that the parameters are achievable in practice, but also that the approximations 
we made in deriving the dispersive Hamiltonian [equation ([8])] remain valid. 

The basic assumptions underlying our description are that the energy scales of the 
TLS and the resonator and the mechanical system are widely separated, hw/A 1 
and that we can only expect to achieve rather weak electro-mechanical coupling, 
K = X/hw <C 1. Furthermore, we assume that the coherence time of the TLS in 
the absence of the resonator, T2, is of order 0.5 /is, in line with recent experimental 
results [22l[23] for a Cooper-pair box embedded in a superconducting cavity. In line with 
this value, we assume a maximum value of ti for the echo experiment of Tc ~ 0.2 fis. 
For concreteness we assume a TLS energy separation 2A/h = 5 GHz and a mechanical 
frequency u!/2tt = 50 MHz. 

Within the regime where uj <^ A/h the maximum amplitude of the mechanical 
motion for which the dispersive Hamiltonian remains valid is set by the condition 
{Xx/xzpAy <^ 1, which we can express as 5 = {2K\ao\hjuj / A)^ <^ 1. We note in passing 
that if |ao| is small enough to satisfy this condition then in practice it will also be 
small enough to ensure that non-linear effects are unimportant in the dynamics of the 
mechanical resonator [13] . 

The value of the electro- mechanical coupling constant, k, which can be achieved of 
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course depends on the actual system used in an experiment. For the specific system we 
have considered here consisting of a mechanical resonator formed by suspending part 
of the qubit circuit [H], the beam is assumed to have a width and thickness ~ 200 nm 
and will need to have a length of a few microns in order to have a frequency of 50 MHz. 
For such a beam Xzp ~ 10"^"^ m and hence we estimate [I6] that coupling strengths up 
to K ~ 0.2 should be achievable. 

The phase space separation which is achieved after a time is S{tc) ^ 2|a;o|^iTc 
(neglecting damping of the mechanical resonator). Using the constraint on the 
magnitude of |q;o|, we obtain S'max — 27r5^/^/t(z/„rc). Assuming (somewhat arbitrarily) a 
value of 5 = 0.04, we find that the maximum value of ao that can be achieved without 
violating our assumptions will be ^ 5/k. Thus for = 0.2 /is and k, = 0.2, we find 
that the maximum value of ao is 25 and S'max = 2.5. This value for the phase space 
separation is encouragingly large, as the minimum uncertainty in phase space of an 
oscillator state with n = 10 (which corresponds to a temperature of about 25 mK for a 
mechanical frequency of 50 MHz) is 2.3. 

6. Results 

We now use the results of the previous section to explore the behaviour of the oscillations 
in during an echo experiment using practicable values of all the parameters. We 
start by examining the envelope of the oscillations in P|_,_^ during an echo experiment 
before and after an inversion pulse at t = ti. The envelope of the oscillations is defined 
by 

i^[^+)W] = ^ + V*/^^e-M^W]. (51) 

where 9(t) is given by equation for times t <ti and by equation (1A.36P for t > ti. 
An example of the expected behaviour as a function of t is shown in figure [2l We assume 
throughout the parameter values discussed in the previous section [uj jl-a = 50 MHz, 
Ua = 5 GHz and T2 = 0.5 /zs) and consider the maximal coupling n = 0.2 and amplitude 
ao = 25. The strength of the mechanical dissipation is specified by the resonator's 
Q-factor, Q = uj/'y. We have taken rT = 10 and as well as considering the case where 
m = n, we also (for theoretical interest) consider the extreme case where the resonator 
is somehow pre-cooled to its ground state, m = 0. 

From the curves in figure [2] we can see that the mechanical resonator is likely to 
have a strong effect on the TLS. It is interesting to compare the curves with and without 
the inclusion of a finite Q-factor for the mechanical resonator. In an echo experiment, 
only mechanical dissipation leads to a deviation from the uncoupled value of at 
the echo point, t = 2ti (i.e., the recoherence). Although an initial mixture of resonator 
states leads to an average over phases associated with each of the different states and 
hence a strong enhancement of the apparent dephasing of the TLS during the first part 
of the experiment {t < ti), after the echo each of these phases unwinds and hence they 
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Figure 2. Envelope of oscillations in in an echo experiment with a tt pulse 
applied at t{— ti) = 0.2 /xs. The blue curves are for k = 0.2, n = rn = 10 and ao = 25. 
The red curves are for the same parameters but with m = 0. In each case the full 
curve is for Q = 3000 and the dashed curve is for the case without any mechanical 
dissipation. The black curve is the result that would be obtained without any coupling 
to the mechanical resonator. 

do not affect tlie behaviour at t = 2ti. On the other hand, when dissipation is included 
we see that the echo signal can be substantially reduced. 

It is important to note that dissipation of the mechanical resonator has only a 
very small effect on the behaviour of the signal £'[P|+)] before the vr pulse is applied. 
This is because the decay of this signal is dominated by the separation of the resonator 
states and the averaging over the different phases associated with each of the states 
in the thermal mixture. The decoherence of the mechanical resonator only starts to 
occur once a superposition has been produced and by the time it has started to develop, 
the value of i?[P|+)] is already close to 0.5. Thus, the decoherence of the mechanical 
resonator can only really be measured by using the echo signal around t ~ 2ti. 

It is interesting to note that pre-cooling the resonator does not affect the echo 
signal by very much. This is again because the phase averaging that occurs for a 
mixed state is largely removed by the use of the echo sequence. However, in the 
presence of dissipation the states involved in a thermal mixture will have slightly different 
amplitudes (compared to the average ao) and hence will all be affected slightly differently 
by the coupling to the environment during the evolution: the mixed initial state curve 
(m = n) does not exactly match the pre-cooled (pure) one (m = 0) at t = 2ti. This 
behaviour can be seen more clearly in figure E] which focusses on the echo signal ait = 2ti 
for a range of values. Over the relatively short time of the echo ti -C I/7 = Q/ur, 
energy diffusion is a very weak effect and hence the evolution of the thermal state is very 
similar to an average over pure initial states with a range of ao values (~ m^/^). Thus 
the results for the initially mixed (m = n) and pure states (m = 0) become very close 
for larger ao values where the variation of the envelope signal with ao is approximately 
linear (on a scale ~ m^/^), and overall the curves are closer for lower n. 
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Figure 3. Envelope of oscillations in in an echo experiment measured at time 
tf — 2ti as a function of ao- The full (dashed) curves are for m = n (to — 0) with 
K = 0.2, ti = 0.2 fis and Q = 10^. The red curves are for n = 20 and the blue curves 
are for n — 10. 




Figure 4. Envelope of oscillations in P\^^ in an echo experiment with a tt pulse applied 
at t{= ti) = 0.2 fis, measured at time tf — 2ti. The blue curves are for k — 0.2, with 
cko = 10 and n = m varied from to 25. The red curves are for the same parameters 
but with n = TO = 10 and ao varied from to 25. In each case the dashed curve is for 
Q — 10"^ and the full curve is for Q = W^. The black line is the result that would be 
obtained without any dissipation to the mechanical resonator. 

In figure HI we compare the effects of varying tlie temperature of tlie meclianical 
resonator's environment and tlie amplitude of tlie initial state on the echo signal at 
t = 2ti. Increasing the value of either n or ao reduces the recoherence at the echo, but 
the dependences are rather different. An important part of any experiment would be 
to test this behaviour, something which could readily be done for ao by simply varying 
the initial drive applied to the mechanical resonator to prepare it in states of different 
amplitude. 



Probing the quantum coherence of a nanomechanical resonator: Echo scheme 19 




Figure 5. Envelope of oscillations in in an echo experiment with a tt pulse 
applied at i = ii, measured at time tf — 2ti as a function oi tf. The full (dashed) 
curves are for k — 0.1 (k = 0.2), with ao — 25 and n = m = 10. The red curves are for 
Q = 10'^ and the blue curves are for Q — lO"'. The black line is the result that would 
be obtained without any dissipation to the mechanical resonator. 

Finally, in figure O we explore how changing the time between the pulses ti (and 
hence the total time for the echo experiment tf = 2ti) affects the behaviour at the 
echo point. This plot shows clearly the strong deviation from simple exponential decay 
that the coupling to the resonator can lead to. As we have already discussed, the 
superposition of resonator states takes time to develop and hence it takes a while before 
decoherence of the mechanical resonator can start to affect the dynamics of the TLS 
which is measured; all the curves in figure [5] initially lie very close to each other. However, 
at longer times the dissipative effect of the mechanical system's environment starts to 
have an important influence. Furthermore, it is clear that for strong enough coupling 
the decay of £'[i^+)(t = 2ti)] occurs on a much faster scale than the relaxation of the 
resonator's energy, 'j = uj/Q, a clear sign that it is the loss of the mechanical system's 
quantum coherence which drives the process. 

The range of Q factors which we have used here, 10'^ — lO'', is appropriate for a 
resonator formed by a suspended metal film |44j . However, where the resonator consists 
of a semiconductor beam which is then coated in a metal layer, somewhat higher Q- 
factors can occur [25] (up to ~ 10^). For very high Q-factor resonators the amplitude 
of the echo signal will be completely dominated by the qubit decoherence and the 
contribution from the resonator's bath may eventually become too small to measure 
in practice. In this regime the measurement of the qubit recoherences would only allow 
an upper bound for the decoherence of the mechanical system to be established. 

7. Conclusions and Discussion 

In this paper we have discussed how a superconducting qubit can be used to probe 
the quantum coherence of a nanomechanical resonator using methods very similar to 
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those applied in recent optical cQED experiments. In particular, we explored how an 
echo experiment could be used to systematically explore the quantum dynamics of a 
mechanical resonator using a superconducting qubit tuned to the degeneracy point as a 
probe. 

The advantages of the echo approach go beyond the practicalities of the system. 
The ability to control the duration of the experiment and to vary the separation of 
resonator states produced (by varying the initial amplitude ao) will make it much easier 
to draw strong conclusions about the nature of the mechanical system's environment. 
Interestingly, we found that over a range of temperatures (corresponding to thermal 
occupation numbers of the resonator up to ~ 20) the recoherences were likely to be 
affected only very weakly by the variance of the initial resonator state implying that it 
is by no means necessary to prepare the resonator in a pure state to obtain important 
information about its quantum dynamics. We expect the echo technique to be rather 
robust in the sense that it should give useful information about the quantum coherence of 
the resonator for a rather wide range of parameters. The larger the separation of states 
achieved during an echo experiment, the more the magnitude of the recoherences will 
tell us about the coherence properties of the mechanical system. However, there is no 
threshold below which nothing is learnt: even if only a very small separation (S < 1) is 
achieved then some information is nevertheless obtained about the dissipative dynamics 
of the mechanical resonator beyond just the energy relaxation rate. 

Since a great deal will be inferred from the deviations between the measured 
dynamics and the reversible dynamics calculated using the dispersive Hamiltonian, it 
will in practice be necessary to be able to discriminate between contributions arising 
from the resonator's environment and those due to the inevitable corrections to the 
model Hamiltonian which is an approximate form. Therefore, an important future 
extension of the current work would be to carry out a systematic numerical study of the 
coupled qubit-resonator dynamics using the full Hamiltonian of the system. Such an 
approach would not just allow us to calculate the effects of corrections to the dispersive 
Hamiltonian, but also allow a more comprehensive modelling of the qubit's environment 
to include energy relaxation. As recent experiments [23] have begun to approach the 
regime where T2 > Ti, the inclusion of a finite Ti is becoming increasingly relevant. 
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Appendix A. Calculation of TLS decoherence for a damped resonator 

In this Appendix we calculate the dynamics of the Wigner function component W-^ 

including the effects of the environment. We start from the equation of motion for 
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[equation fl36p ]. which in terms of the interaction picture, 

p+_(t)=e[^2A/m/T.]t^_^_(^)^ 

becomes 

= —i \uja)a, P+-] — iuoi {a^a + 1/2, p+--} + 
Defining the Wigner transform in the usual way, 

W+.{x,p]t) = ^ 

we obtain the Wigner-transformed equation of motion. 



dy{x + y\p+^{t)\x-y)e-''^y'^ 
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x' 
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where x = x/x^p and p = p{2/mhioY^'^. 

In order to solve this equation of motion we make a Gaussian ansatz, assuming 
that the Wigner function takes the form of a Gaussian multiplied by a phase factor 



W+-{x,p;t) = WG{x,p;t)— = — -= 



(A.5) 



where D is the determinant of the matrix 

^xp ^p 



(A.6) 



and the five parameters (x, p, ax-, CTp, (Jxp) and the phase 6 are taken to be time dependent. 
Defined in this way Wg{x,p) is normalized (i.e. integrating it over all x,p values gives 
unity) and so the factor of has been introduced as Tr[p+_(t)] is by definition 

(for a TLS) a complex number with amplitude < 1/2. The initial Gaussian remains a 
Gaussian for all times (albeit with different parameters) and hence remains normalized, 
thus 



Tr[p+_(t)] = dp dS:WG{S:,p;t) 



and hence 



Tr[p+_(t)] =e-^2AtM-VT.! 



(A.7) 



(A.8) 



This function is all we need to calculate the probability of finding the qubit in state |+), 

1 



;i-2Im{Tr[p+_(t)]}) 



(A.9) 



Thus using the definition of the initial state of the TLS [equation [9] we can see that 
^'(0) = 37r/2. 

In principal we can solve for the time dependence of the six parameters in the 
Gaussian by substituting the ansatz into the equation of motion directly [57] and 



(A.4) 
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equating powers of x,p. However, in practice the problem is more readily solved 
using the characteristic function which is defined by the relation |l5] 



G{k,q) 



dx / dpW+-{x,p;t)e 



ikxiqp 

e 



2 



(A.IO) 



(A.ll) 



The equation of motion for the characteristic function is readily derived from the 
corresponding one for the Wigner function, 



dG , , dG , , dG 



(A.12) 



[7(n + l/2)+zc.i/2] {e + q^)G + ^ 



d'^G d^G 

+ 



dk"^ dq^ 

Substituting the trial function into the left-hand side of this equation and equating 
powers of /c, g, kq etc, leads directly to a set of equations of motion for the six time 
dependent parameters. 
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(A.13) 

(A.14) 
(A.15) 
(A.16) 

xp ^. (A. 17) 

,p - 7[(Tp- AT] -z^i((Tj + (T^,p- 1), (A.18^ 

where we have defined A^ = 2?T + 1. We now need to solve these equations subject to 
appropriate initial conditions. 

Assuming a thermal state displaced by the coherent amplitude ao, the set of initial 
conditions is as follows: x(0) = 2Re[ao]! p(0) = 2Im[ao], cr^(O) = crp(O) = 2m + 1 
and o"^p(0) = 0. With these initial conditions it is clear that Gxp will remain zero 
for all times and the position and momentum variances will always remain the same, 
= crp(^) = 1 + oit), following the simplified equation 

a = — 7[cr — cTg] — iuJi{2a + cr^) 

where (Tq = cr(0) = 2m. The solution of this equation gives 
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The final part of the calculation involves calculating x{t) and p{t) and hence 
obtaining the phase 0{t). The equations for the averages are most easily solved in 
terms of the variables ai(2) = {x + {—)ip)/2 which obey the equations of motion 



a2= yiuj iui{l + cr(t)) j 02 

and can be integrated to give, 

ai(t) = ai(0)e(-^--'^/2)tg-i./o'[i+<x{t0]rft' 
a2(t) = a2(0)e(*"-^/2)tg-.c./*[i+a(O]dt'. 

The integral in the exponentials is readily evaluated. 



[1 + a{t')]dt' 
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2uJi 



and hence we find 

ai(2)(t) = ai(2)(0)e(-(+)^™^^)* 
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(A.22) 

(A.23) 
(A.24) 

(A.25) 
(A.26) 



The initial values of ai(2) are ai(0) = ao, 02(0) = q;q. 

Finally then we are in a position to obtain the required phase, 0'{t). Noting that 
-|-p2 _ 4q^q2 and using the appropriate initial condition (^'(0) = 37r/2), we obtain 



e'{t) = e'{0) -ioi [ (1 + a{t'))dt' - 2uJi [ ai{t')a2{t')dt' 
Jo Jo 

= 3n/2 



t^ + uji/3)t 



— i- 



;i-M) 



i In 



-i2ujil3t 



1-M 



I _ ]\/fQ-i2uj^pt 
- 1 



(A.27) 
(A.28) 



Thus we arrive at our final result 
Tr[p+-(t)] = (y 



-i2At/h-t/T2+i9(t) 



(A.29) 



where we have defined 6{t) = 6'(t) — 6'{0). This result [equation lA.29] and the expression 
for 6'{t) above gives equation fj42|) in the main text. 

We now extend this calculation to consider the spin-echo case where the system is 
prepared and allowed to evolve in the way we have been considering, but after time ti 
an additional control pulse is applied to invert the populations of the two eigenstates. 
The system is then allowed to evolve for a further time t2 before a final control pulse is 
applied and then a measurement is made. 

In order to obtain p+-(t/ = ti + ^2) we need to solve equation (1A.2P twice: first 
for the period ti and then using the Hermitian conjugate of this solution as the initial 
condition for a further evolution over time ^2- As before we use the Wigner function 
approach and hence use [W^+-(ii)]* as an initial condition for equation flA.4p . 
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Using the above calculation we can immediately write down 



-WG{x,p]ti), 



where the Gaussian Wigner function is in this case parameterized by 
(yx{ti) = o-p(ii) = o"i(^i) + 1 



(Tl(ti 



7 



2a;i 

ai(ti) = a*e*("+"^^*)*i 



I _|_ ]\^*gj2aji/3*ti 
I _ _/\^*gi2a;i/3*ti 

1 - M* 



]_ _ ]\/[*Qi2uJi/3*ti 
1-M* 



with 



37r/2 - ( y ) ti -iln 



1-M* 



X _ ]\/[*Qi2uJi(3*ti 



(A.30) 

(A.31) 
(A.32) 

(A.33) 

(A.34) 

(A.35) 



— t- 



|ao| 



/3* 



'1-M*] 



\ _ ]\^*Qi2uJif3*ti 

The final step is then to use W^__{ti) as the initial condition for evolved over a 
time ^2- Solving the equations of motion for the Gaussian parameters [equations flA.lQp . 
( ]A.2ip and ( ]A.22p ] using the initial conditions given by equations ( ]A.32p -( lA.34[) above 
we finally obtain the phase parameter which is used in equation (!5T!) for t > ti, 

1-M' 



e{tj) = + 37r/2) -['^ + u:il3]t2-i\n 



1^ 



Q-i2uiiPt2 _ ]^ 
X _ ]\/['Q-i2uJil3t2 



where 



M' 



1 



(A.36) 



(A.37) 
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